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We study the occurrence and nature of naked singularities for a dust model with non-zero cos- 
mological constant in (n + 2)-dimensional Szekeres space-times (which possess no Killing vectors) 
for n > 2. We find that central shell-focusing singularities may be locally naked in higher dimen- 
^ ' sions but depend sensitively on the choice of initial data. In fact, the nature of the initial density 

determines the possibility of naked singularity in space-times with more than five dimensions. The 
^-ji results are similar to the collapse in spherically symmetric Tolman-Bondi-Lemaitre space-times. 

Ctf ■ 

PACS numbers: 0420D, 0420J, 0470B 

o : 
<n : 

! I. INTRODUCTION 

> ; 

An extensive study [1-6] of gravitational collapse has been carried out of Tolman-Bondi- 
Lemaitre (TBL) spherically symmetric space-times containing irrotational dust. Due to 
, simplifications introduced by the spherical symmetry several generalizations of this model 

have been considered. A general conclusion from these studies is that a central curvature 
singularity forms but its local or global visibility depends on the initial data. Also the study 
of higher-dimensional spherical collapse reveals the interesting feature that visibility of 
\ singularity is impossible in space-times with more than five dimensions with proper choice 

O -1 °f regular initial data [7-9] . 

*-i ■ 

By contrast, there is very little progress in studying non-spherical collapse. The basic 
difficulty is the ambiguity of horizon formation in non-spherical geometries and the influence 
• i-h , of gravitational radiation. Schoen and Yau [10] proposed a sufficient criterion for the 

■ formation of trapped surfaces in an arbitrary space-time but it fails to say anything about 

the conditions which lead to the formation of naked singularities. This problem has been 
restated by Thorne [11] in the form of a conjecture known as the hoop conjecture which 
states that "horizon form when and only when a gravitational mass M gets compacted into 
a region whose circumference in every direction is C < AirM". Subsequently, there were 
attempts namely numerical simulations of prolate and oblate collapse [12], gravitational 
radiation emission in aspherical collapse [13], analytical studies of prolate collapsing 
spheroids [14] and others [15,16] to prove or disprove the conjecture. Interestingly all of 
them either confirmed or failed to refute the conjecture. 

The quasi-spherical dust collapse models, given by Szekeres metric [17] were analyzed 
by Szekeres himself [18], Joshi and Krolak [19], Deshingkar, Joshi and Jhingan [20] and 
extensively by Goncalves [21]. The solutions for dust and a non-zero cosmological constant 
were found by Barrow and Stein-Schabes [26]. In this work, we study the gravitational 
collapse in the recently generalized (n + 2)-dimensional Szekeres metric [22]. As in 
four-dimensional space-time, this higher dimensional model does not admit any Killing 
vector and the description quasi-spherical arises because it has invariant family of spherical 
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hypersurfaces. The paper is organized as follows: In section II we derive the basic equations 
and regularity conditions. In section III we investigate the formation and the nature of 
central singularity. We study the formation of an apparent horizon due to collapse in 
section IV. The nature and the strength of the singularity is investigated by an analysis 
of the geodesies in sections V and VI respectively. Finally the paper ends with a short 
discussion. 



II. BASIC EQUATIONS AND REGULARITY CONDITIONS 

Recently, dust solutions have been obtained for (n + 2)-dimensional Szekeres' space-time 
metric for which the line element is [22] 



ds 2 = dt 2 - e 2a dr 2 - e 2 ^ dx i W 

»=i 

where a and j3 are functions of all the (n + 2) space-time coordinates. Under the assumption 
that (= -g^) 7^ 0, the explicit form of the metric coefficients are [22] 

e p = R{t,r) e »(r,xi,..,x n ) ( 2 ) 

and 



R' + R v' 



(3) 



where 



e-^4(r)^^+^B 1 (r) I ,+C(r) (4) 



i=i i=i 
and R satisfied the differential equation 



Here A is cosmological constant, /(r) and F(r) are arbitrary functions of r alone; and 
the other arbitrary functions, namely A(r), Bi(r) and C(r), are restricted by the algebraic 
relation [22] 



Y^B 2 -AAC = -1 (6) 



The r-dependence of these arbitrary functions A, Bi and C plays an important role in 
characterizing the geometry of the (n + l)-dimensional space. In particular, if we choose 



A(r) = C(r) = 5 and Bi{r) = (V i = 1,2, ...,n) then using the transformation [22] 

xi = Sin6 n Sin9 n -x Sin92Cotj9\ 

X2 = Cos6 n Sin6 n -i Sin92Cot\9\ 

x 3 = Cosd n -iSin9 n -2 Sin02Cot\9\ 



x n -\ — Cos9 3 Sin9 2 Cot^9i 
x n = Cos8 2 Cot^9i 
the space-time metric (1) reduces to the spherically symmetric TBL form: 

ds 2 = dt 2 - - R '\ dr 2 - R 2 dn 2 n . (7) 
l + /(r) 

In the subsequent discussion we shall restrict ourselves to the quasi-spherical space-time 
which is characterized by the r dependence of the function v (i.e., v' ^ 0). From the Einstein 
field equations we have an expression for energy density for the dust, as 

N n F' + (n + l)iV 
P (t,r, Xl ,...,x n ) = - Rn{RI + Rul) • (8) 

Thus a singularity will occur when either (i) R = i.e., (3 = — oo or (ii) a = — oo. Using 
the standard terminology for spherical collapse, the first case corresponds to a shell-focusing 
singularity, while the second case gives rise to a shell-crossing singularity. As in a TBL 
space-time, the shell-crossing singularities are gravitationally weak, and we shall not con- 
sider them any further here. Hence, we shall restrict ourselves to the situation with a > — oo. 

Suppose that f = t, is the initial hypersurface from which the collapse develops. For the 
initial data we assume that R(ti,r) is a monotonically increasing function of r. So, without 
any loss of generality, we can label the dust shells by the choice R(ti,r) — r. As a result, 
the expression for the initial density distribution is given by 

x n F' + (n+l)Fv' 
Pi(r,xi,...,x n ) = p(U,r,x 1 ,...,x„) = - r n( 1 + n/ ;) ( 9 ) 

If we started the collapse from a regular initial hypersurface the function p i must be 
non-singular (and also positive for a physically realistic model). Furthermore, in order for 
the space-time to be locally flat near r — 0, we must have f(r) — > as r — > 0. Then, from 
equation (5), the boundedness of R 2 as r — > demands that F(r) <~ 0(r m ) where m>n—l. 
But, for small r, p^r) ~ ^iL an d consequently, for regular p^r) near r — 0, we must have 
F(r) ~ 0(r n+1 ). Thus, starting with a regular initial hypersurface, we can express F(r) 
and Pi (r) by power series near r = as 



F(r)=J2Fjr n+j+1 (10) 

3=0 



and 



3=0 



(11) 
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As v 1 appears in the expression for the density as well as in the metric coefficient, we can 
write 



i=-i 



(12) 



where u_ 1 > —1. 



Now, using these series expansions in equation (9) we have the following relations between 
the coefficients, 



n(n + 1) n ( 1 



P2= 2 



n+ 1 + 



F 2 



P 3 = 2 



n + 1 + — ) F 3 



1 + 



i 7 ! 



(13) 



and so on. 



Finally, if we assume that the density gradient is negative and falls off rapidly to zero near 
the centre then we must have p 1 = and p 2 < 0. Consequently, we have the restrictions 
that Fi = and F 2 < 0. 



III. FORMATION OF SINGULARITY AND ITS NATURE 

In order to form a singularity from the gravitational collapse of dust we first require that 
all portions of the dust cloud are collapsing i.e., R < 0. Let us define t s f(r) and t sc (r) as the 
time for shell-focusing and shell-crossing singularities to occur occurring at radial coordinate 
r. This gives the relations [18] 



R(t sf ,r) = 



and 



R'(t sc ,r) + R(t sc ,r)i/'(r,x 1 ,x 2 , -,x n ) = 0. 
Note that 'i sc ' may also depend on x\, x%, x n . 



(14) 



(15) 



As mentioned earlier, the shell-crossing singularity is not of much physical interest so we 
shall just consider the shell-focusing singularity for the following two cases: 

(i) Marginally bound case : f(r) = with A ^ 

In this case equation (5) can easily be integrated to give 



2n 



(n + 1)A 



Sink 1 



2Ar n+1 



n(n+ l)-F(r) 



Sink' 1 



/ 2KR n + 1 
n(n + l)F{r) 



(16) 



Also, from equations (14) and (16), we have 
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2n 



(n + l)A 



Sink 



/ 2Ar" +1 
n(n + l)F(r) 



(17) 



(ii) Non-marginally bound case with time symmetry : f(r) ^ 0, R(ti,r) = 0, A = 
Here, the solution of equation (5) gives t as a function of r : 



t = n 



r(6+l/2) 2 ll 2' ' ' ' 



{n + l)VF 

and so the time of the shell-focusing singularity is given by 



(18) 



2v^F r(6+l) n+1 

tsf ~ U+ (n + l)^r(6+l/2) r 2 ' 



(19) 



where 2-F1 is the usual hypergeometric function with 6 = | + . However, for the 
five-dimensional space-time (n = 3), i? has the particularly simple form 



R 2 = r*-^{t-uf, 



(20) 



and therefore 



tsf =U + 



(21) 



IV. FORMATION OF TRAPPED SURFACES 

The event horizon of observers at infinity plays an important role in the nature of 
the singularity. However, due to the complexity of the calculation, we shall consider a 
trapped surface which is a compact space-time 2-surface whose normals on both sides 
are future-pointing converging null-geodesic families. In particular, if (r = constant, t = 
constant) the 2-surface 5 r> t is a trapped surface then it, and its entire future development, 
lie behind the event horizon provided the density falls off fast enough at infinity. 

If t a h is the instant of the formation of apparent horizon then we must have [7,8,18] 

2AR n+1 (t ah , r) - n(n + l)^" 1 ^, r) + n(n + l)F(r) = 0. (22) 
Thus, for the above two cases, the explicit expressions for t a h are the following: 



t = U 



2n 



_ (n + l)A 
for f(r) = 0, A ^ and 



Sink 1 



2Ar n+1 



i(n+l)F(r) 



Sink- 1 



l 2AR"+i(t ah ,r) 
n(n + l)F(r) 



(23) 



tah — ti H~ 



F 



T(b+ 1/2) 



(24) 



G 



for /(r) 7^ 0, R(ti,r) = 0, A = and for n — 3 we get 



t a h — U 



(25) 



From the expressions for t s f and t a h, we note that the shell-focusing singularity that 
appears at r > is in the future of the apparent horizon in both cases. But since we 
are interested in the central shell-focusing singularity (at r = 0) we require the time of 
occurrence of central shell-focusing singularity t (= t s f(0)). From equations (16), (19) and 
(22), taking the limit as r — > 0, we have that 



t a h(r) - to 



1 



.'i/2 



2A 



n(n+l)F 



Fr 



F 2 



(4A + 3n(n+l)F )F 2 . 2 
4F (2A + n(n+l)F ) 1 ' 



_ (4A + 3n(n+l)F )F 1 F 2 (32A 2 + 40n(n + 1)AF + 15n 2 (n + l) 2 + Fg)F? \ , 

^3 ; ^ 77. , — ; 7 , i \ t— i \ 1 . — o — , ; tt-^tt; I T 



4F (2A + n(n+ 1)F ) 



24F 2 (2A + n(n + l)F ) 



+. 



(n + l)^" 1 
for /(r) = 0, A ^ and 



n + l 1 

F r^ + p ir - 

n — 1 



+ 



n- 1 



n - 2 Fi \ sn-i 

■ 1 r Ti^r + ... 



2(n - 1) F a 



(26) 



r(6 + i) 



(n+l)F 3/2 r(6+l/2) 



Fir 



4F y 



Fi 



9FjF 2 
4F 



5Fi 3 \ 



(« + W 



F r" 



+ ■ 



n 



-Fir » 



+ 



(27) 



for /(r) 7^ 0, r) = 0, A = and 

for n = 3 (with /(r) ^ 0, i?(t i; r) = 0, A = 0) 



tah{r) - t = -- 



2F 



3/2 



Fx r + 



F2 + f °~Jf 



r 2 + . 



(28) 



Here, to is the time at which the singularity is formed at r = 0, and t a h(r) is the instant 
at which a trapped surface is formed at a distance r. Therefore, if t a h(r) > to, a trapped 
surface will form later than the epoch at which any light signal from the singularity 
can reach an observer. Hence the necessary condition for a naked singularity to form is 
t a h(f) > to; while that for a black hole to form is t a h{r) > to- It is to be noted that this 
criteria for naked singularity is purely local. Hence, in the present problem it possible to 
have local naked singularity or a black hole form under the conditions shown in the Table I. 



TABLE-I 





Naked Singularity 


_ffo/e 




Marginally bound case: 
/(r) = 0,A^0 


(i) p x < 0, Vn 

(ii) Pl = 0,p 2 <0,n = 2 

(iii) p! = 0, p 2 < 0, n = 3, 
F 2 <-2F o yi + 4 

(iv) p! = 0, p 2 = 0, p 3 < 0, 
n = 2, 


(i) p 1 =0,p 2 <0,n = 3, 
F 2 >-2F o yi + 4 

(ii) Pi = 0,p 2 < 0,n> 4 

(iii) Pl = 0, p 2 = 0, p 3 < 0, n = 2, 

i^3>-2Fo i/2 v /l+3| r 

(iv) pj = 0, p 2 = 0, p 3 < 0, n > 3 

(v) pj = 0,p 2 = 0, ..^pj = 0, 
Pi+i < 0, j > 3,Vn 




F 3 <-2^ / yi+^- 


Non-marginally bound 
case with time symmetry: 
/(r)^0, R(ti,r)=0, ' 
A = 


(i) Pi < 0,Vn 

(ii) p x = 0,p 2 < 0,n = 2 

(iii) p! = 0, p 2 < 0, n = 3, 

(iv) p! = 0, p 2 = 0, p 3 < 0, 


(i) p 1 = 0,p 2 <0,n = 3,F 2 >-F 2 

(ii) pj = 0,p 2 < 0,n > 4 

(iii) pj = 0, p 2 = 0, p 3 < 0, n = 2, 

(iv) pj = 0, p 2 = 0, p 3 < 0, n > 3 

(v) pj = 0,p 2 = 0, ...,p j = 0, 
Pj+i < 0,j > 3,Vn 





However, if we relax the condition on the density gradient (namely the condition that 
p x = 0, p 2 < i.e., Fi — 0, F 2 < 0) then it is possible to have a naked singularity in 
all dimensions. These results are identical to those obtained in the presence of spherical 
symmetry. Hence the local nature of singularity is not affected by the geometry of the 
space-time with respect to whether it is spherical or quasi-spherical. However, we note that 
the Szekeres space-times retain certain special geometrical properties shared by spherical 
space-times in all dimensions, in particular they do not allow gravitational radiation to be 
present in the space-time. 

V. GEODESICS AND THE NATURE OF SINGULARITY 

Here, for simplicity, we shall consider only the marginally bound case (i.e., f(r) = 0) with 
A = 0. Then R(t,r) has the simple explicit solution (choosing the initial time t% = 0): 

2 

(29) 

We now follow the geodesic analysis of Joshi and Dwivedi [23] for TBL model. For 
simplicity of calculation we introduce the following functions: 



R 



n+i n + 1 rpn~\ , 
' 2 2~V F (r) t 
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where the constant a is restricted by a > 1. In order to determine the nature of the 
singularity we examine whether it is possible to have outgoing null geodesies which are 
terminated in the past at the central singularity r = 0. Suppose this occurs at singularity 
time t = to at which R(t , 0) = 0. In order to decide this we start with radial null geodesic, 
given by 



dt _ a _ R'Q - RQ' 



Next, we introduce the notation u = r a , so we write 



^ = U(X,u), (32) 



where U(X, u) = (l - y/^j f + l^g^ with H = ^ X + . 

We now study the limiting behaviour of the function X as we approach the singularity at 
R = 0, u = along the radial null geodesic identified above. If we denote the limiting value 
by Xq then 



X = lim f = Urn f§ = lira U(X,u) =U(X ,0) 

R^Ou^O R^Ou^O R^Ou^O 



(33) 



Thus if this polynomial equation in X has at least one positive real root then it is possible 
to have a radial null geodesic outgoing from the central singularity. More explicitly, X is 
the root of the following equation in X: 



( a_ ^Tl) X " + (^TT _% ) ^ XH± -0oX^ + v /^ e = (34) 

where suffix 'o' for the variables £,?y, ( and 6 stands for the values of these quantities at 
r = 0. As exact analytic solution is not possible for X so we study the roots by numerical 
methods. Table II shows the dependence of the nature of the roots on the variation of the 
parameters involved. 

Table II: Positive roots (Xo) of the eqn. (34) for different values of the parameters 
namely a, rj 07 ^ Q7 ( , <d and in different dimensions, n. 
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a 


Vo 


£o 


Co 


e 






Positive roots (X ) 
















4D 


5D 


6D 


7D 


8D 


10D 


14D 


1 


-6 


.05 


.01 


-5 


.00997 


.09864 


.21131 


.30963 


.38979 


.50797 


.64774 


2 


2 


5 


3 


2 


3.3019, 
3 


1.732, 
1.633 


1.4422, 
i .oiyo 


1.1968, 

1 31 fi1 
1 .0101 


1.2457, 

1 1 3/L^ 
1 .1040 


1.1699, 
1.075 


1.105, 
1.0334 


3 


2 


5 


3 


2 


- 


- 








- 


- 


5 


4 


.1 


3 


2 


1.7068, 
.77462 


1.2678, 
.80678 


1 . 1010, 

.82742 


111 A9. 

1 . 1140, 

.84873 


i .uooy, 
.86583 


1.0612, 
.89102 


1.0377, 
.92113 


1 
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1 


1 


2.035 


1.6116 


I.OIOO 


l.OUZo 


l.OOll 






1 





4 


1 


1 


14.8239 


- 








- 


- 


1 


.1 


1 


1 


1 






- 


- 


- 






1 


1 


1 


1 


1 


1.3104, 
1 


1.1547, 
1 


1.0934, 

i 
i 


1.0627, 
i 


1.045, 

i 
i 


1.0265, 
1 


1.0124, 
1 


4 





5 


.1 


4 
















1 





1 


1 


9 


2.25, 
1.5963 


2.5336, 
1.1725 


2.2120, 
1.1045 


1.9675, 
1.0752 


1.7967, 
1.0588 


1.5829, 
1.041 


1.3756, 
1.0256 


1 





1 


5 


1 
















1 





1 


.1 


1 


.90297, 
.10733 


.87007, 
.36226 


.86041, 
.54185 


.85671, 
.66031 


.85235, 
.74455 







From the definition of the functions in equation (30) it is clear that the parameters a, £ , £ 
are always positive while r] , Q can have positive as well as negative values. From the Table 
we see that as we decrease the value of a, keeping the other parameters fixed, then the 
formation of a naked singularity is more probable in higher dimensions. The situation for £ 
is similar. But for £ and rj the condition is reversed: black hole formation becomes more 
probable as we increase the dimension of the space-time if we decrease the value of £ (or 
77 ) and keep the other constants fixed. However, no such conclusions can be drawn about 
the variation of 6q. 



VI. STRENGTH OF THE NAKED SINGULARITY 

A singularity is called gravitationally strong, or simply strong, if it destroys by crushing 
or tidally stretching to zero volume all objects that fall into it; it is called weak if no object 
that falls into the singularity is destroyed in this way. A precise characterization of Tiplcr 
strong [24] singularities has been given by Clarke and Krolak [25] , who proposed the strong 
focusing condition. A sufficient condition for a strong curvature singularity is that, for at 
least one non space-like geodesic with affine parameter A — > on approach to the singularity, 
we must have 

Urn \ 2 R l:j K l K3 > ,„-,. 
A -> (6b) 

where K 1 = is the tangent vector to the radial null geodesic. 
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Our purpose here to investigate the above condition along future-directed radial null 
geodesies that emanate from the naked singularity. Now equation (35) can be expressed as 
(using L'Hopital's rule) 



lira XiRijKW ^(H^x^-in+i)^) 

' ' ' Co 



A^O 2X^N +r, 

where H = H(X ,0),N = N(X Q ,0). 

The singularity is gravitationally strong in the sense of Tipler if 

(n+l)6 l 



£ - (n + l)r/ Q > max < 0, -- 



n+l 

x 2 



or 



^ - (n + l)r} < min \ 0, - — + }}®° 



x 



(36) 



If the above condition is not satisfied for the values of the parameters then 
^ _^ q X 2 RijK t K : ' < and the singularity may or may not be Tipler strong. 



VII. DISCUSSION AND CONCLUDING REMARKS 

In this paper we have studied gravitational collapse in (n + 2)-dimensional space-time 
using a higher-dimensional generalization of the quasi-spherical Szekeres metrics with non- 
zero cosmological constant. We have examined the local nature of the central shell-focusing 
singularity by a comparative study of the time of the formation of trapped surface and the 
time of formation of the central shell-focusing singularity. If we assume the initial density 
gradient falls off rapidly and vanishes at r = 0, then naked singularity formation is possible 
only up to space-time dimension five. However, if we drop the above restriction on the 
initial density distribution then the Cosmic Censorship Conjecture may be violated in any 
dimension (n > 2). Thus we deduce that the nature of the central singularity depends 
sensitively in these metrics on the choice of the initial data (particularly on the choice of 
initial density profile). This is also confirmed by our geodesic study following the approach 
of Joshi and Dwivedi [23] where we have shown numerically that the nature of singularity 
depends on the value of the defining parameters at r = 0. Finally, we have examined 
the strength of the naked singularity using the criterion introduced by Tipler [24]. We 
found that the naked singularity will be a strong curvature singularity depending on the 
appropriate choice of the value of the parameters at r = 0. These investigations provide 
some insights into the phenomenon of gravitational collapse in a situation without any 
imposed Killing symmetries. However, the collapses are special in other senses which permit 
exact solutions to be found. In particular, there is an absence of gravitational radiation in 
these space-times [27]. An investigation of its role is a challenge for future analytic and 
computational investigations. 
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